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ABSTRACT 

The  characteristics  and  dynamics  of  inertia-gravity  waves  generated  in  the  vicinity  of  an  intense  jet  stream/ 
upper-level  frontal  system  on  18  February  2001  are  investigated  using  observations  from  the  NOAA  Gulfstream- 
IV  research  aircraft  and  numerical  simulations.  Aircraft  dropsonde  observations  and  numerical  simulations 
elucidate  the  detailed  mesoscale  structure  of  this  system,  including  its  associated  inertia-gravity  waves  and  clear- 
air  turbulence.  Results  from  a  multiply  nested  numerical  model  show  inertia-gravity  wave  development  above 
the  developing  jet/front  system.  These  inertia-gravity  waves  propagate  through  the  highly  sheared  flow  above 
the  jet  stream,  perturb  the  background  wind  shear  and  stability,  and  create  bands  of  reduced  and  increased 
Richardson  numbers.  These  bands  of  reduced  Richardson  numbers  are  regions  of  likely  Kelvin-Helmholtz 
instability  and  a  possible  source  of  the  clear-air  turbulence  that  was  observed. 


1.  Introduction 

Gravity  waves  are  ubiquitous  in  the  atmosphere  and 
represent  a  signihcant  process  within  the  momentum 
budget  of  the  middle  atmosphere  (see,  e.g.,  Fritts  and 
Alexander  2003).  Gravity  waves  can  be  generated  by 
mountains,  convection,  and  baroclinic  waves.  Compo¬ 
nent  structures  that  exist  within  baroclinic  waves,  such 
as  jets  and  upper-level  fronts,  are  important  sources  of 
gravity  waves  and  lower-frequency  inertia-gravity 
waves.  Upper-level  fronts  have  also  been  identihed  as 
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a  preferred  location  for  clear-air  turbulence  (CAT)  (e.g., 
Kennedy  and  Shapiro  1980).  CAT  remains  a  serious 
concern  for  the  aviation  industry,  and  can  also  affect 
the  mixing  and  stratosphere-troposphere  exchange  of 
chemical  constituents  (e.g.,  Shapiro  1980;  Ravelin  et  al. 
2002).  The  relationship  between  upper-level  fronts  and 
CAT  has  mostly  been  attributed  to  Kelvin-Helmholtz 
instabilities  related  to  frontogenesis  or  flow  deformation 
(e.g.,  Keller  1990).  Also,  frontally  forced  gravity  waves 
have  been  identified  as  a  possible  source  of  CAT  (e.g., 
Mancuso  and  Endlich  1966;  Knox  1997);  however,  lim¬ 
ited  progress  has  been  made  in  validating  this  hypoth¬ 
esis.  This  study  investigates  the  generation  of  inertia- 
gravity  waves  by  a  jet  stream/upper-level  front,  and  the 
causal  relationship  between  the  inertia-gravity  waves 
and  CAT. 

Previous  studies  have  investigated  the  generation  of 
inertia-gravity  waves  within  idealized  fronts  (e.g.,  Sny- 
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der  et  al.  1993;  Griffiths  and  Reeder  1996),  simple  bar- 
oclinic  waves  (e.g.,  O’Sullivan  and  Dunkerton  1995; 
Zhang  2004),  or  through  observations  near  jets  (e.g.. 
Guest  et  al.  2000;  Plougonven  et  al.  2003).  Also,  there 
have  been  a  number  of  real  case  studies  (e.g..  Powers 
and  Reed  1993;  Zhang  et  al.  2001)  focusing  on  the  role 
of  tropospheric  inertia-gravity  waves  in  contributing  to 
the  organization  of  convection  and  affecting  weather. 
The  waves  in  the  aforementioned  studies  possessed  a 
variety  of  scales,  ranging  from  horizontal  wavelengths 
of  0(100  km)  (mesoscale  inertia-gravity  waves),  up¬ 
ward  of  0(1000  km). 

Inertia-gravity  waves  have  been  observed/simulated 
in  all  regions  of  baroclinic  waves,  but  most  of  the  pre¬ 
vious  studies  have  identihed  regions  within  the  baro¬ 
clinic  wave  that  were  preferable  for  the  formation  of 
inertia-gravity  waves.  For  example,  using  observations 
Guest  et  al.  (2000)  showed  that  the  region  between  an 
upstream  trough  and  downstream  ridge  is  a  likely  region 
for  inertia-gravity  waves  in  the  lower  stratosphere.  This 
region  was  also  the  focus  of  Zhang’s  (2004)  simulations 
of  an  idealized  baroclinic  wave.  O’Sullivan  and  Dunk- 
erton’s  (1995)  simulation,  however,  highlighted  inertia- 
gravity  waves  generated  in  the  jet  streak  how  between 
an  upstream  ridge  and  downstream  trough.  Nonetheless, 
although  these  preferred  regions  of  wave  activity  have 
been  identihed,  the  existence  of  waves  in  other  parts  of 
the  baroclinic  wave  is  not  unlikely.  The  case  considered 
herein  is  in  the  northwesterly  how  of  a  jet  streak,  in  the 
Northern  Hemisphere,  between  an  upstream  ridge  and 
downstream  trough. 

The  mechanisms  responsible  for  the  excitation  of  these 
inertia-gravity  waves  remains  a  subject  of  research.  The 
proposed  mechanisms  include  adjustment  from  an  un¬ 
balanced  to  a  balanced  state  (geostrophy  or  a  higher-order 
balance)  (e.g.,  Uccellini  and  Koch  1987;  Zhang  2004), 
and  nonlinear  forcing  analogous  to  Lighthill  radiation 
(Reeder  and  Griffiths  1996).  All  of  these  proposed  ex¬ 
citation  mechanisms  characteristically  take  place  in  re¬ 
gions  of  strong  jets  associated  with  baroclinic  waves. 

The  specihc  details  of  the  above  wave  source  mech¬ 
anism  are  not  the  focus  of  this  study.  Instead,  this  study 
focuses  on  a  jet  stream/upper-level  front  that  was  ob¬ 
served  during  the  Severe  Clear-Air  Turbulence  Collid¬ 
ing  with  Aircraft  Traffic  (SCATCAT)  experiment. 
SCATCAT  was  an  experiment  designed  to  examine  CAT 
associated  with  jet  streams  and  upper-level  fronts.  The 
National  Oceanographic  and  Atmospheric  Administra¬ 
tion’s  (NOAA)  Gulfstream-IV  (G-IV)  weather  recon¬ 
naissance  aircraft  was  the  principal  observing  platform, 
encompassing  high-resolution  dropsonde  releases  and 
in  situ  measurements.  On  the  17-18  February  2001 
hight,  the  G-IV  observed  an  intense  jet  stream  and  an 
upper-level  front,  and  encountered  moderate  CAT.  The 
turbulence  was  encountered  in  a  region  of  possible  grav¬ 
ity  wave  activity. 

This  paper  has  two  objectives.  The  hrst  objective  is 
to  describe  the  detailed  mesoscale  structure  of  the  jet/ 


front  system  that  was  observed  surrounding  0000  UTC 
18  February  2001.  The  second  objective  is  to  utilize 
hnescale  numerical  simulations  to  examine  the  gener¬ 
ation  of  inertia-gravity  waves  by  the  jet/front  system 
and  relate  these  waves  to  small-scale  instabilities,  to 
advance  our  knowledge  of  the  role  of  inertia-gravity 
waves  in  generating  CAT.  It  will  be  shown  that  the 
simulations  compare  qualitatively  well  to  the  observa¬ 
tions,  especially  considering  the  paucity  of  data  in  the 
Pacihc  and  the  scales  of  motion  involved.  Although  the 
simulations  possess  some  limitations,  they  provide  use¬ 
ful  insight  into  the  dynamical  behavior  of  the  observed 
system  and  the  relationship  between  inertia-gravity 
waves  and  CAT. 

The  remainder  of  the  paper  is  organized  as  follows. 
The  observations  from  the  G-IV  surrounding  0000  UTC 
18  February  2001  are  presented  (section  2),  followed 
by  a  discussion  of  the  synoptic-scale  how  derived  from 
a  forecast  (section  3).  In  section  4,  a  nested  model  is 
used  to  examine  the  hnescale  structure  of  the  how,  in¬ 
cluding  the  jet,  upper-level  front,  and  inertia-gravity 
waves.  It  is  shown  that  these  waves  are  responsible  for 
generating  a  Kelvin-Helmholtz  instability.  In  section  5, 
the  NOAA  G-IV  dropsondes  are  analyzed  for  inertia- 
gravity  waves  to  determine  whether  the  inertia-gravity 
waves  derived  from  the  simulations  are  realistic.  Finally, 
the  results  are  summarized  in  section  6. 

2.  Dropsonde  observations:  17-18  February  2001 

Surrounding  0000  UTC  18  February  2001,  the  G-IV 
observed  an  intense  jet  stream/frontal  system  over  the 
Pacihc  Ocean  at  approximately  40°N,  200°E.  During 
this  hight,  17  closely  spaced  dropsondes  were  released 
between  2326  UTC  17  February  2001  and  0024  UTC 
1 8  February  200 1  with  an  average  horizontal  separation 
of  approximately  40  km.  The  hight  track  of  the  G-IV 
during  the  dropsonde  deployment  is  shown  in  Fig.  1; 
the  G-IV  was  initially  hying  to  the  north,  then  turned 
to  the  northeast.  Also  shown  in  Fig.  1  is  the  approximate 
location  of  the  jet  stream  at  0000  UTC  18  February, 
derived  from  a  24-h  forecast  from  the  Naval  Research 
Laboratory’s  Coupled  Ocean-Atmosphere  Mesoscale 
Prediction  System  (COAMPS;  see  Hodur  1997). 
(CO AMPS  and  its  conhguration  are  described  in  section 

3. )  The  17  dropsonde  prohles  were  combined  to  form 
a  vertical  cross-section  analysis.  The  close  separation 
of  these  sondes  along  with  their  high  vertical  resolution 
permits  the  analysis  of  coherent  structures.  In  order  to 
construct  this  cross-section  analysis,  the  prohle  from 
each  sonde  was  interpolated  onto  a  regular  vertical  grid 
with  50-m  spacing.  These  data  were  then  hltered  to 
remove  signals  with  vertical  scales  of  motion  less  than 
250  m.  Each  sonde  was  assigned  a  distance  along  the 
hight  track  (relative  to  the  hrst  sonde  released)  using 
its  release  position.  The  resulting  cross  sections  of  po¬ 
tential  temperature  and  wind  speed  are  shown  in  Eig. 
2.  Three  of  the  17  sondes  did  not  record  winds.  At  the 
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Fig.  1.  The  G-IV  flight  track  during  the  dropsonde  deployment 
(thick  line)  and  the  release  point  of  every  second  dropsonde  (crosses). 
Position  of  the  jet  stream  at  0000  UTC  18  Feb  2001  from  a  COAMPS 
24-h  forecast,  illustrated  by  the  9  km  AMSL  wind  vectors  and  wind 
speed  contoured  and  shaded  at  50,  60,  and  70  m  s”'.  Also  shown  is 
the  location  of  COAMPS  domain  2,  and  the  CHM  domain  1. 

top  of  each  sounding,  there  is  one  data  point  at  the 
release  height  determined  from  the  G-IV  in  situ  flight 
data,  and  because  it  takes  a  finite  time  for  the  sonde  to 
reach  equilibrium  with  its  environment,  there  is  a  layer 
of  missing  data  about  1  km  thick  below  flight  level.  The 
data  are  interpolated  across  this  layer  of  missing  data. 

Figure  2  shows  the  detailed  structure  of  the  jet/front 
system.  The  observed  maximum  wind  speed  of  the  jet 
exceeds  100  m  s  ',  with  its  core  at  approximately  9  km 
above  mean  sea  level  (AMSL).  There  is  also  evidence 
of  a  weaker  Jet  (65  m  s“')  at  approximately  5.5  km 
AMSL.  The  frontal  zone  below  the  jet  is  between  ap¬ 
proximately  2  and  6  km  AMSL,  and  tilts  at  about  100:1 
(horizontakvertical).  Above  the  jet  are  perturbations  in 
potential  temperature  that  may  be  due  to  gravity  waves. 
These  perturbations  are  in  the  region  just  to  the  cyclonic 
(northeast)  side  of  the  jet  core,  extending  to  the  anti- 
cyclonic  side  (southwest).  A  line  that  represents  a  pos¬ 
sible  tilting  phase  structure  is  also  shown  overlaid  in 
Fig.  2a.  Above  and  below  the  jet  core,  the  vertical  shear 
is  large  and  exceeds  approximately  20  m  s^'  km“‘  (0.02 
s“').  The  negative  (speed)  shear  above  the  jet  is  gen¬ 
erally  stronger  than  the  positive  (speed)  shear  below. 

After  deploying  the  dropsondes,  the  G-IV  completed 
stacked  flight  legs  in  the  region  of  possible  gravity 
waves  above  the  jet  core.  During  these  stacked  legs  the 
G-IV  encountered  moderate  CAT.  Regions  of  CAT  that 
induced  vertical  accelerations  of  the  G-IV  that  exceeded 
4.9  m  s“^  (0.5  g’s)  in  magnitude  are  also  shown  in  Fig. 
2a.  These  regions  of  turbulence  coincide  with  the  strong 
negative  shear  region  above  the  jet,  and  the  perturba¬ 
tions  in  potential  temperature  that  may  be  gravity  wave 
induced.  A  detailed  examination  of  the  turbulence  and 
in  situ  data  from  the  G-IV  can  be  found  in  Koch  et  al. 
(2003)  and  Koch  et  al.  (2004,  manuscript  submitted  to 
J.  Atmos.  Sci.). 


(a)  Potential  Temperoture  (K) 


(b)  Wind  Speed  (m  s’') 


Fig.  2.  Cross  sections  of  (a)  potential  temperature  and  (b)  wind 
speed  derived  from  the  dropsondes.  The  horizontal  axes  run  from  the 
southwest  to  the  northeast  along  the  flight  track  shown  in  Fig.  1.  The 
locations  of  the  dropsondes  with  good  thermodynamic  and  wind  data 
are  shown  by  the  downward-pointing  aiTows.  The  upward-pointing 
arrow  marks  the  change  in  direction  of  the  G-IV  shown  in  Fig.  1. 
The  thick  dashed  lines  mark  the  flight  level  and  first  level  of  good 
data.  Also  shown  in  (a)  is  a  line  that  represents  a  possible  gravity 
wave  phase  line,  and  regions  where  the  G-IV  encountered  vertical 
accelerations  that  exceeded  4.9  m  s"^  in  magnitude  (diamonds). 


3.  The  synoptic-scale  flow 

COAMPS  is  used  to  analyze  and  forecast  the  evolving 
three-dimensional  flow  surrounding  the  observations. 
This  numerical  model  is  a  finite-difference  approxi¬ 
mation  to  the  fully  compressible,  nonhydrostatic  equa¬ 
tions  and  uses  a  terrain-following  vertical  coordinate 
transformation.  The  finite-difference  schemes  are  of  sec¬ 
ond-order  accuracy  in  time  and  space.  Reflection  of 
waves  at  the  upper  boundary  is  mitigated  by  a  radiation 
condition  following  Klemp  and  Durran  (1983)  and  Bou- 
geault  (1983).  The  simulations  use  a  full  suite  of  phys¬ 
ical  parameterizations  to  represent  boundary  layer,  ra¬ 
diative,  and  moist  processes  (e.g.,  Hodur  1997;  Hodur 
and  Doyle  1999). 

For  this  experiment,  COAMPS  is  configured  with  two 
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(a)  Surface  Pressure  and  9km  Wind  Speed 


Fig.  3.  (a)  Surface  pressure  (mb)  and  (b)  the  800-600-mb  thickness 
(dam)  at  0000  UTC  18  Feb  2001,  from  the  COAMPS  (domain  1) 
24-h  forecast.  Also  shown  in  both  plots  is  the  9  km  AMSL  wind 
speed  contoured  and  shaded  at  50,  60,  and  70  m  s“'. 


(a)  Potential  Temperature  (K)  0000  UTC  18  Feb  2001 


Fig.  4.  Cross  sections  of  (a)  potential  temperature  and  (b)  wind 
speed  from  COAMPS  (domain  2)  at  0000  UTC  18  Feb  2001.  The 
cross  section  is  between  the  southwest  and  northeast  corners  of  the 
model  domain  (see  Fig.  1).  Also  shown  in  (a)  is  the  2-PVU  contour 
of  potential  vorticity  (thick  line). 


horizontally  nested  domains  of  97  X  73  and  133  X  133 
points,  with  horizontal  grid  spacings  of  54  and  18  km, 
respectively.  The  model  upper  boundary  is  at  30  km 
with  80  vertical  levels,  with  a  vertical  grid  spacing  of 
10  m  at  the  lowest  level  stretched  to  400  m  aloft.  The 
forecast  is  initialized  at  0000  UTC  17  February  2001. 
The  initial  fields  for  the  model  are  created  from  mul¬ 
tivariate  optimum  interpolation  analyses  of  upper-air 
sounding,  surface,  aircraft,  and  satellite  data  that  are 
quality  controlled  and  blended  with  the  12-h  CO  AMPS 
forecast  fields.  Lateral  boundary  conditions  make  use 
of  the  Navy  Operational  Global  Analysis  and  Prediction 
System  forecast  fields  following  Davies  (1976). 

Figures  1  and  3  show  the  synoptic  features  from  do¬ 
main  1  of  CO  AMPS  at  0000  UTC  18  February  2001 
(24-h  forecast).  At  this  time,  the  jet  that  is  the  focus  of 
this  study  flows  from  the  northwest,  downstream  of  the 
surface  high  and  upstream  of  the  surface  cyclone  (Fig. 
3a).  An  upper-level  front  is  associated  with  this  jet.  This 
front  is  aligned  with  the  jet,  as  shown  by  the  800-600- 
mb  thickness  (Fig.  3b).  As  expected  from  the  thermal 


wind  balance,  the  front  has  strong  gradients  in  potential 
temperature  perpendicular  to  the  jet  on  its  cyclonic 
(northeast)  side. 

The  vertical  structure  of  the  jet/front  system  is  shown 
in  southwest  to  northeast  cross  sections  at  0000  and 
0600  UTC  (Figs.  4  and  5,  respectively).  At  0000  UTC 
(Fig.  4)  the  jet  maximum  is  at  about  8  km  AMSL  and 
the  maximum  wind  speed  exceeds  80  m  s  '.  There  are 
strong  horizontal  shears  on  the  cyclonic  side  of  the  jet, 
and  the  isopleths  are  nearly  vertical.  Below  the  jet  is 
the  upper-level  front  with  its  strong  horizontal  and  ver¬ 
tical  gradients  in  both  wind  speed  and  potential  tem¬ 
perature.  Also,  the  height  of  the  tropopause  descends 
significantly  from  its  background  height  of  about  10- 
14  km,  into  a  deep  fold  on  the  cyclonic  side  of  the  jet, 
following  the  front  to  about  4  km.  On  the  anticyclonic 
side  of  the  jet  core,  the  static  stability  in  the  upper 
troposphere  is  very  low  over  a  depth  of  about  3  km  and 
500-1000  km  in  length. 

At  0600  UTC  (Fig.  5),  the  jet/front  system  has  a  similar 
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(a)  Potential  Temperoture  (K)  0600  UTC  18  Feb  2001 


Structure  to  that  of  0000  UTC  (Fig.  4).  However,  the  jet 
is  faster  and,  subsequently,  the  vertical  shears  above  and 
below  and  the  horizontal  shear  on  the  cyclonic  side  of 
the  jet  are  greater.  The  frontal  baroclinicity  intensifies 
and  steepens.  The  tropopause  fold  maintains  its  vertical 
extent  but  has  broadened  in  the  horizontal.  Of  special 
importance  to  this  study  are  strong  perturbations  in  the 
potential  temperature  that  extend  vertically  into  the 
stratosphere  above  the  jet  core.  These  perturbations  are 
the  signature  of  vertically  propagating  inertia-gravity 
waves  with  horizontal  wavelengths  (in  this  southwest  to 
northeast  plane)  of  about  250  km.  These  waves  are  dis¬ 
cussed  in  further  detail  in  section  4. 

The  simulated  structure  of  the  jet/front  system  (Figs. 
4  and  5)  is  qualitatively  similar  to  that  shown  in  the 
dropsonde  analysis  (Fig.  2).  However,  there  are  some 
differences;  for  example,  the  speed  of  the  jet  is  under¬ 
represented,  and  the  simulated  upper  front  develops 
about  6  h  later  than  observed.  Also,  the  modeled  upper 
front  is  steeper  than  the  observations  and  tilts  at  about 
60;  1  (horizontakvertical).  These  differences  are  ex¬ 
pected  considering  the  model’s  location  in  the  data- 
sparse  Pacific  and  are  not  a  major  concern.  In  the  model 
the  amplitude  of  the  inertia-gravity  waves  at  0600  UTC 


is  largest.  Therefore,  even  though  our  observations  sur¬ 
round  0000  UTC,  we  have  selected  0600  UTC  as  the 
focus  of  our  further  modeling  efforts  and  analyses,  and 
our  comparisons  to  the  observations  will  be  qualitative. 

4.  Finescale  modeling 

In  the  previous  section,  the  synoptic  forecast  pro¬ 
duced  inertia-gravity  waves  above  the  jet/upper-front 
system.  In  this  section,  a  finescale  nested  model  (the 
Clark-Hall  model)  is  used  to  better  resolve  the  upper 
front,  inertia-gravity  waves,  and  possible  CAT  gener¬ 
ation. 

a.  The  Clark-Hall  model  and  nesting  procedures 

The  Clark-Hall  model  (CHM;  Clark  1977;  Clark  et 
al.  1996)  is  a  nonhydrostatic,  anelastic  finite  difference 
mesoscale  model.  In  this  study,  it  is  configured  in  three 
dimensions,  warm-rain  cloud  processes  are  treated  using 
the  Kessler  (1969)  scheme,  and  the  subgrid  mixing  is 
the  Smagorinsky  (1963)-Lilly  (1962)  closure.  Radiative 
and  surface  fluxes  are  neglected.  The  CHM  is  capable 
of  two-way  nesting  (Clark  and  Hall  1991),  which  in¬ 
cludes  refinement  in  both  the  horizontal  and  vertical,  as 
well  as  at  elevated  domains.  This  feature  of  the  model 
allows  small-scale  features  to  be  examined  in  significant 
detail,  while  still  allowing  larger-scale  features  to  evolve 
realistically. 

One  other  feature  of  the  CHM  is  its  ability  to  be  (one 
way)  nested  within  an  operational-style  model.  In  this 
study,  COAMPS  provides  initial  conditions  to  the  CHM. 
At  subsequent  times,  the  COAMPS  model  fields  are 
used  as  boundary  conditions  for  the  CHM  outermost 
domain.  These  boundary  conditions  are  constructed  us¬ 
ing  linear  interpolation  from  hourly  COAMPS  data.  [A 
more  detailed  description  of  the  way  the  CHM  is  forced 
by  a  larger-scale  model  is  found  in  Clark  et  al.  (2000).] 

The  horizontal  locations  of  the  CHM  domains  used 
in  this  study  are  shown  in  Fig.  6,  and  their  geometries 
are  summarized  in  Table  1 .  (Note,  each  model  grid  uses 
constant  Alatitude  and  Alongitude  spacing,  and  there¬ 
fore  although  AT  is  constant,  A  A  is  a  function  of  lati¬ 
tude.  AA  =  AT  cos(latitude)/cos(41.5),  where  41.5°N 
is  the  latitude  at  the  center  of  the  outermost  domain.) 
Not  all  domains  are  used  in  all  experiments  and  we 
denote  experiment  n  as  a  model  simulation  that  uses  n 
domains  only.  The  vertical  grid  spacing  in  domain  1 
was  chosen  to  match  that  of  CO  AMPS;  this  configu¬ 
ration  reduces  mismatches  between  COAMPS  and  the 
CHM  near  the  lateral  boundaries  of  CHM  domain  1. 
The  uppermost  10  km  of  CHM  domain  1  includes  a 
Rayleigh  friction  absorber  to  reduce  reflection  of  dis¬ 
turbances  from  the  upper  boundary.  Domain  3  has  the 
same  horizontal  grid  spacing  as  domain  2,  but  smaller 
vertical  grid  spacing.  This  domain  configuration  pro¬ 
vides  a  smoother  transition  from  domain  1  to  3  by  re¬ 
ducing  inner-domain  boundary  effects,  and  allowing  an 
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Fig.  6.  The  horizontal  location  of  CHM  domains  2-4  within  domain 
1  (see  Fig.  1).  Also  shown  is  a  grid  (dashed);  each  element  of  the 
grid  has  dimensions  of  100  km  X  100  km. 

effective  4;  1  refinement  in  the  vertical.  Domain  3  is  used 
for  most  of  the  analysis  of  the  jet/front  system,  and 
domain  4  is  used  to  examine  shearing  instabilities  in  the 
region  of  the  inertia-gravity  waves.  Although  there  are 
only  a  few  clouds  in  the  CHM  domain,  the  explicit 
microphysics  of  the  CHM  (i.e.,  no  cumulus  parameter¬ 
ization)  restricts  the  horizontal  grid  spacing  of  domain 
1 .  The  horizontal  grid  spacing  of  6  km  was  found  to  be 
about  the  largest  grid  spacing  that  did  not  produce  prob¬ 
lematic  results  due  to  improper  treatment  of  cloud  evo¬ 
lution. 

b.  High-resolution  simulation  with  the  Clark-Hall 
model 

Figure  7  presents  a  cross  section  of  simulated  potential 
temperature  (Fig.  7a)  and  wind  speed  (Fig.  7b)  from 
CHM  domain  3  at  0600  UTC  for  experiment  3.  The 
potential  temperature  identihes  the  front,  which  is  stron¬ 
gest  between  4  and  8  km.  The  strongest  part  of  the  front 
is  associated  with  a  deep  tropopause  fold  that  extends 
downward  to  4  km.  The  background  tropopause  height 
ranges  from  5  km  on  the  northeast  side  of  the  front  to 
10  km  on  the  southwest  side.  In  addition  to  the  deep 
tropopause  fold,  there  is  a  shallow  secondary  fold  (at 
around  400-km  distance  and  8  km  AMSL),  and  undu¬ 
lations  in  the  tropopause  on  the  anticyclonic  side  of  the 
secondary  fold.  These  undulations  have  a  horizontal  scale 
of  around  150  km.  The  wind  speed  shows  the  jet  core 
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Fig.  7.  Cross  sections  of  (a)  potential  temperature  and  (b)  wind 
speed  from  CHM  (domain  3,  experiment  3)  at  0600  UTC  18  Feb 
2001.  The  cross  section  is  between  the  southwest  and  northeast  cor¬ 
ners  of  the  model  domain  (see  Fig.  6).  Also  shown  in  (a)  is  the  2- 
PVU  contour  of  potential  vorticity  (thick  line). 

to  be  at  about  9  km  AMSL,  and  through  this  cross  section 
the  maximum  wind  speed  is  about  75  m  s  ‘  (although 
the  maximum  wind  speed  in  the  domain  is  about  85 
m  s  ').  The  vertical  wind  shear  is  maximized  above  the 
jet  and  within  the  frontal  zone  associated  with  the  tro¬ 
popause  fold.  In  these  regions,  the  vertical  shear  exceeds 
20  m  s^’  km^'  (0.02  s^').  In  addition  to  these  vertical 
shears,  the  jet  exhibits  two  regions  of  localized  cyclonic 
shear  at  the  level  of  maximum  wind  (9  km)  at  distances 
of  ~425  and  —550  km.  These  horizontal  shear  zones  are 
within  the  secondary  and  deep  tropopause  folds.  The 


Table  1.  Geometry  of  the  CHM  domains:  the  meridional  grid  spacing  (AT),  vertical  grid  spacing  (AZ),  lateral  dimensions  (NX  =  NY), 
vertical  dimensions  (NZ),  lowermost  domain  level  (Z,,),  uppermost  domain  level  (Z„),  and  domain  initialization  time  T,,. 


Domain 

AT  (km) 

AZ  (m) 

NX 

NZ 

Z„  (km) 

Z„  (km) 

To 

1 

6 

400 

222 

92 

0 

33.6 

1800  UTC  17  Feb 

2 

3 

200 

222 

120 

0 

21.2 

0000  UTC  18  Feb 

3 

3 

100 

180 

160 

0 

13.3 

0000  UTC  18  Feb 

4 

1 

50 

200 

100 

8 

12.9 

0400  UTC  18  Feb 
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(a)  Potential  Temp,  (K)  0300  UTC  18  Feb  2001 


(b)  Potential  Temp.  (K)  0600  UTC  18  Feb  2001 


X(km) 


Fig.  8.  Horizontal  cross  sections  of  potential  temperature  through 
6  km  AMSL  at  (a)  0300  and  (b)  0600  UTC  from  CHM  (domain  3, 
experiment  3).  The  contour  interval  is  0.5  K. 


overall  shape  of  the  jet  is  similar  to  that  analyzed  from 
the  dropsonde  observations;  however,  the  modeled  jet  is 
weaker  and  therefore  the  modeled  vertical  shears  above 
and  below  the  jet  are  also  weaker  than  observed. 

The  simulation  shows  evidence  of  vertically  propa¬ 
gating  inertia-gravity  waves,  which  appear  as  wavelike 
perturbations  in  the  potential  temperature  above  the  jet 
core,  extending  to  its  anticyclonic  side  (see  Fig.  7a).  Like 
the  COAMPS  forecast,  these  waves  appear  after  0000 
UTC  and  their  amplitude  increases  a  few  hours  before 
0600  UTC;  their  amplitudes  are  maximized  at  0600  UTC. 


Potential  Temp.  (K)  0600  UTC  18  Feb  2001 


Fig.  9.  Horizontal  cross  section  of  potential  temperature  through 
1 1  km  AMSL  at  0600  UTC  from  CHM  (domain  3,  experiment  3). 
To  assist  with  interpretation,  the  potential  temperature  is  contoured 
as  the  perturbation  from  335  K  with  negative  values  dashed.  The 
contour  interval  is  0.5  K.  Also  shown  are  lines  that  approximate  lines 
of  constant  gravity  wave  phase. 


This  increase  in  wave  amplitude  coincides  with  rapid 
upper-level  frontogenesis.  The  intensification  of  the  front 
is  illustrated  in  Fig.  8,  which  shows  the  6-km  potential 
temperature.  At  0300  UTC,  the  upper  front  is  forming, 
with  the  strongest  horizontal  gradients  at  its  northwestern 
end.  In  time,  these  gradients  progressively  increase  along 
the  front  from  the  northwest,  until  0600  UTC  when  the 
front  has  intensified  with  strong  horizontal  gradients 
along  its  entire  length.  The  mature  front  is  oriented  par¬ 
allel  to  the  jet  from  the  northwest  to  southeast. 


c.  Analysis  of  inertia-gravity  waves  from  the 

simulation 

The  inertia-gravity  waves  in  Fig.  7a  appear  as  a  quasi- 
monochromatic  wave  packet  with  approximately  two 
horizontal  wavelengths  evident.  The  phase  lines  of  this 
wave  tilt  toward  the  southwest,  implying  vertical  prop¬ 
agation.  Using  Fig.  7a,  the  vertical  wavelength  of  this 
wave,  in  the  lower  stratosphere,  is  estimated  to  be  ap¬ 
proximately  2.5  km.  Therefore,  Iml  =  2.513  X  10“^ 
rad  m  ',  and  assuming  upward  energy  propagation  and 
positive  intrinsic  frequency  m  is  negative. 

Figure  9  illustrates  the  horizontal  structure  of  the 
waves  using  the  11 -km  AMSL  potential  temperature 
from  domain  3  of  experiment  3.  Also  shown  in  this 
figure  are  two  lines  that  represent  lines  of  constant  in¬ 
ertia-gravity  wave  phase,  passing  through  adjacent 
wave  maxima  and  minima.  The  amplitude  of  the  wave- 
induced  perturbations  in  potential  temperature  is  ap- 
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proximately  3  K.  Using  Fig.  9  the  zonal  wavelength  is 
estimated  to  be  approximately  240  km,  and  the  merid¬ 
ional  wavelength  is  approximately  280  km.  The  hori¬ 
zontal  wave  vector,  =  {k,  /),  is  perpendicular  to  the 
lines  of  constant  phase,  and  the  direction  of  the  vertical 
tilt  shown  in  Fig.  7a  implies  that  the  horizontal  wave 
vector  points  to  the  southwest.  Therefore,  Ky  = 
(-2.618,  -2.244)  X  10-"  rad  m-‘,  IKhI  =  3.448  X 
10“^  rad  m  ‘,  and  the  horizontal  wavelength  Ay  ~  180 
km.  This  estimate  of  is  consistent  with  measurements 
from  Figs.  7a  and  9. 

The  dispersion  relation  for  nonhydrostatic  gravity 
waves  with  rotation  is 

+  F  ~  0)2  -  /"’ 


ambiguously  whether  the  mesoscale  structures  are  the 
source  of  the  gravity  waves  or  a  response  to  some  other 
forcing  that  is  also  generating  the  waves.  Such  source 
analysis  is  beyond  the  scope  of  this  study,  but  is  an  area 
of  active  research. 

Einaudi  and  Lalas  (1975)  showed  that  a  wave  prop¬ 
agating  through  an  almost  neutral  atmosphere  can  sig- 
nihcantly  affect  the  local  Richardson  number  (Ri)  and 
possibly  induce  shearing  instabilities.  In  the  case  pre¬ 
sented  here,  the  atmosphere  is  stably  stratihed,  yet  in  the 
region  above  the  jet  core  the  wind  shear  is  sufficiently 
strong  that  the  Richardson  number  is  uniformly  low  (Ri 
~  2).  The  Richardson  number  is  dehned  as  N^/ 1  Def  1 2, 
where  N  is  the  Brunt-Vaisala  frequency  and  I  Def  I  is 
the  magnitude  of  the  total  deformation,  but  at  the  scales 
of  motion  considered  herein. 


where  N  is  the  Brunt-Vaisala  frequency,  w  is  the  in¬ 
trinsic  frequency  of  the  gravity  wave,  and  /  is  the  Cor¬ 
iolis  parameter.  Making  the  hydrostatic  approximation 
(^2  +  I-  m2),  the  dispersion  relation  can  be  written  as 


Therefore,  using  the  estimates  of  k,  I,  and  m  made  earlier, 
N  =  0.02  rad  s  '  (a  value  representative  of  the  modeled 
lower  stratosphere)  and  /  =  10  rad  s  ',  the  intrinsic 
frequency  is  approximately  2.952  X  10  ''  rad  s  '  and 
the  intrinsic  period  is  approximately  5.9  h.  These  wave 
characteristics  are  extremely  similar  to  those  in  other 
modeling  studies  such  as  those  of  Zhang  (2004);  how¬ 
ever,  the  wave  phase  lines  are  oriented  parallel  to  the 
jet,  whereas  the  waves  that  were  the  focus  of  Zhang’s 
study  were  perpendicular  to  the  jet.  Yet,  F.  Zhang  (2004, 
personal  communication)  notes  that  gravity  waves 
aligned  with  the  northwesterly  flow,  like  those  consid¬ 
ered  here,  are  also  present  in  his  idealized  simulation. 

As  mentioned  earlier,  the  amplification  of  these  in¬ 
ertia-gravity  waves  coincides  with  the  intensification  of 
the  upper  front.  Also,  the  wave  phase  lines  are  approx¬ 
imately  aligned  with  the  upper  front  (e.g..  Figs.  8  and 
9).  These  two  features  of  the  inertia-gravity  waves  im¬ 
ply  that  the  waves  are  indeed  generated  by  upper-level 
frontogenesis.  However,  the  mesoscale  features  of  the 
jet/front  system,  including  the  horizontal  shear  layers 
adjacent  to  the  jet  and  the  associated  deep  and  secondary 
tropopause  folds,  possess  horizontal  scales  similar  to 
the  generated  inertia-gravity  waves.  Without  a  sophis¬ 
ticated  source  analysis  it  is  difficult  to  determine  un- 


where  u  is  the  zonal  velocity  and  v  is  the  meridional 
velocity,  is  a  good  approximation.  If  the  horizontal  ve¬ 
locity  components  and  the  potential  temperature  are  sep¬ 
arated  into  a  background  profile  and  a  gravity  wave- 
induced  perturbation,  the  Richardson  number  can  be 
written  as 

gfdl  ^ 

6\dz  dz ) 

Ri  =  ^ ^ - , 

(dU  didV  (dV  St/V 

\  t/z  dz  )  \  t/z  dz  ) 

where  the  overbar  denotes  the  background  quantity  and 
the  prime  denotes  the  perturbations.  Assume  that  the 
wave  perturbations  in  zonal  velocity,  meridional  veloc¬ 
ity,  and  potential  temperature  take  the  form 

u'  =  A„  cos((/)  +  aj),  v'  =  A„  cos(<p  +  Wj),  and 
6'  =  Ag  sm(4>). 

The  wave  phase  is  cj)  =  kx  +  ly  +  mz  —  wt,  and 
are  phase  differences;  and  A„,  A„,  and  Ag  are  the  wave 
amplitudes  in  zonal  velocity,  meridional  velocity,  and 
potential  temperature,  respectively.  [These  amplitudes 
and  phase  differences  are  not  independent  and  are  cou¬ 
pled  via  the  polarization  relations  for  an  inertia-gravity 
wave  (Gill  1982).]  Invoking  the  Wentzel-Kramers-Bril- 
louin  (WKB)  approximation  (essentially  m  =  constant), 
the  Richardson  number  can  be  written  as 


0  dz 


niAg  cos((j)) 


dU 

2 

dV 

- mA„  sin((p  +  aj) 

dz 

-f 

- mAy  sm((p  +  a-,) 

dz 
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Fig.  10.  Richardson  number  (shaded)  at  0600  UTC  18  Feb  2001 
from  CHM  (domain  3,  experiment  3).  The  cross  section  is  from  the 
southwest  to  northeast  corners  of  the  domain  (see  Fig.  6).  Also  shown 
is  potential  temperature  at  2-K  intervals. 

Inertia-gravity  waves  will  induce  variations  in  the  local 
Brunt-Vaisala  frequency  and  vertical  wind  shear,  which 
will  in  turn  affect  the  Richardson  number.  Therefore, 
the  Richardson  number  may  be  signihcantly  affected  by 
gravity-wave-induced  perturbations  if 

.  de  ,  dU  ^  dV 

mA„ - ,  mA„ - ,  or  mA„ - . 

"  dz  "  dz  dz 

In  the  center  of  domain  3  of  experiment  3  at  0600  UTC, 
b^ween  10  and  12  km  AMSL,  dOldz  ~  0.01  K  m^', 
dUldz  ~  —0.01  s^',  and  dV/dz  ~  0.01  s“'.  Also,  from 
the  simulation,  the  amplitudes  of  the  inertia-gravity 
waves  at  11  km  AMSL  are  estimated  as  A„  ~  5  m  s“', 
~  5  m  s  ',  and  A^  ~  3  K.  Therefore,  using  these 
quantities  with  the  vertical  wavenumber  determined  ear¬ 
lier  (m  =  —2.5  X  10  rad  m  '),  all  three  of  the  above 
conditions  are  met,  implying  that  the  inertia-gravity 
waves  may  signihcantly  affect  the  local  Richardson 
number. 

Consider  a  southwest  to  northeast  cross  section  of  Ri 
(=  A^/IDefU)  from  domain  3  (experiment  3)  of  the 
CHM  (Fig.  10).  This  hgure  shows  that  in  the  lower 
stratosphere,  in  the  region  of  vertically  propagating 
waves,  the  Richardson  number  is  perturbed  signihcant¬ 
ly.  Subsequently,  the  Richardson  number  held  possesses 
coherent  bands  of  reduced  Richardson  numbers  (Ri  < 

1)  and  bands  of  increased  Richardson  numbers  (Ri  > 

2) .  These  bands  of  reduced  and  increased  Richardson 
numbers  are  aligned  with  the  gravity  wave  phase  lines. 
The  bands  of  reduced  Richardson  numbers  approxi¬ 
mately  coincide  with  the  wave-induced  reduction  in  sta¬ 
bility,  that  is,  in  quadrature  with  the  potential  temper¬ 
ature  perturbations.  Thus,  as  the  gravity  waves  propa¬ 
gate  vertically  through  the  strong  negative  shear  region 
above  the  jet,  they  perturb  the  background  how  signif¬ 
icantly  and  cause  coherent  bands  of  low  Richardson 
numbers  (<1  and  <1/4  in  some  small  regions).  These 
bands  are  regions  of  possible  Kelvin-Helmholtz  insta¬ 


bility  and  CAT.  Thus,  although  these  waves  are  not  nec¬ 
essarily  strong  in  amplitude,  their  high  vertical  wave- 
number  coupled  with  the  large  background  wind  shear 
and  low  background  Richardson  number  allow  them  to 
induce  a  shearing  instability.  In  an  observational  case 
study.  Ravelin  et  al.  (2001)  also  showed  that  inertia- 
gravity  waves  were  responsible  for  generating  a  shear¬ 
ing  instability.  However,  in  Ravelin  et  al.’s  study  the 
environment  possessed  low  wind  shear  (high  back¬ 
ground  Richardson  number),  but  the  amplitude  of  the 
inertia-gravity  waves  was  large  enough  to  induce  strong 
perturbations  in  wind  shear  and  stability,  reduce  the 
Richardson  number,  and  induce  a  shearing  instability. 

Figure  10  also  shows  a  large  region  of  Ri  <  1  in  the 
upper  troposphere  (altitudes  between  6  and  9  km  AMSL, 
and  between  0-  and  400-km  distance).  This  region  is 
due  to  the  strong  shear  below  the  jet  core  and  the  rel¬ 
atively  low  upper-tropospheric  static  stability.  Other  re¬ 
gions  of  low  Richardson  numbers  also  exist  within  the 
upper  front. 

Figure  1 1  shows  the  simulated  helds  from  domain  4 
at  0600  UTC  18  February  2001,  for  experiment  4,  in 
the  region  above  the  jet  where  the  inertia-gravity  waves 
are  affecting  the  Richardson  number.  The  subgrid  tur¬ 
bulent  kinetic  energy  (TKE)  (Fig.  11a)  is  inferred  from 
the  model  eddy  viscosity,  K^,  using  Deardorff’s  (1980) 
formulation;  TKE''^  =  \0KJI,  where  /  =  AZ  =  50  m. 
Nonzero  subgrid  TKE  exists  in  the  strong  shear  above 
and  below  the  jet,  and  in  the  region  of  the  vertically 
propagating  inertia-gravity  waves.  Like  the  bands  of 
reduced  Richardson  numbers,  the  subgrid  TKE  is  lo¬ 
cated  in  the  region  of  the  wave-induced  reduction  in 
stability,  and  also  in  regions  where  the  potential  tem¬ 
perature  perturbation  is  negative.  The  subgrid  TKE  is 
not  in  exact  quadrature  with  the  perturbation  potential 
temperature  because  of  the  contribution  of  the  (back¬ 
ground  and  wave  perturbed)  wind  shear  to  the  Rich¬ 
ardson  number.  In  this  highest  resolution  domain,  the 
subgrid  TKE  is  not  as  coherent  as  may  be  expected  from 
Fig.  10  because  parts  of  the  region  of  shear  instability 
are  mixed  out  by  the  resolved  and  parameterized  flow. 
Figure  lib  shows  an  estimate  of  the  resolved  TKE 
(RTKE),  where  RTKE  =  («""  -f  v"^  +  w"^)l2.  The 
perturbation  velocities  u",  v",  w"  are  dehned  as  devia¬ 
tions  from  a  smoothed  version  of  the  background  ve¬ 
locities.  (This  smoothing  removes  scales  of  motion  that 
are  less  than  20  km  in  the  horizontal.)  The  resolved 
TKE  held  identihes  the  existence  of  small-scale  re¬ 
solved  motions  (<20  km)  that  are  stronger  than  the 
subgrid  TKE  and  that  are  mixing  out  the  shear  insta¬ 
bilities  in  some  regions.  These  small-scale  motions  are 
evident  as  small  perturbations  in  the  wind  held  (Eig. 
11c). 

The  subgrid  TKE  is  enhanced  in  two  main  regions. 
First,  in  horizontally  coherent  layers  above  and  below 
the  jet  core  in  the  regions  of  strongest  vertical  wind 
shear.  These  turbulent  layers  may  not  be  related  to  the 
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Fig.  11.  Cross  sections  of  (a)  subgrid  TKE,  (b)  resolved  TKE,  and 
(c)  wind  speed  from  CHM  (domain  4,  experiment  4)  at  0600  UTC 
18  Feb  2001.  The  cross  section  is  from  the  southwest  to  northeast 
corners  of  the  domain  (see  Fig.  6).  In  (a)  and  (b)  the  shading  is  linear 
with  the  maximum  values  of  0.2  and  0.7  m^  s”^  (respectively)  marked 
by  darkest  shading.  Wind  speed  contours  in  (c)  are  0.5  m  s“‘.  Also 
shown  in  (a)  and  (b)  is  potential  temperature  at  2-K  intervals. 


Fig.  12.  Horizontal  cross  section  of  subgrid  TKE  through  11 
km  AMSL  from  CHM  (domain  4,  experiment  4)  at  0600  UTC  18 
Feb  2001.  Shading  is  the  same  as  in  Fig.  11a,  also  shown  is  the 
0.01  m^  s”-  contour. 

inertia-gravity  waves,  and  are  not  featured  in  the  re¬ 
solved  TKE  field.  The  second  region  is  within  the  low- 
stability  phase  of  the  inertia-gravity  waves,  where  both 
subgrid  and  resolved  TKE  are  enhanced.  The  resolved 
TKE  may  be  identifying  growing  modes  associated  with 
a  Kelvin-Helmholtz  instability.  However,  despite  our 
efforts  with  nesting,  the  grid  spacing  in  this  domain  is 
still  too  coarse  to  properly  resolve  the  Kelvin-Helm¬ 
holtz  waves  that  are  expected  from  such  a  shear  insta¬ 
bility.  Also,  the  simulation  time  of  this  domain  may  be 
insufficient  to  properly  spin  up  such  motions. 

At  upper  levels  the  subgrid  TKE  shown  in  Eig.  11a 
is  located  along  the  phase  lines  of  the  gravity  waves, 
but  the  structures  are  incoherent  or  patchy  in  the  vertical. 
However,  a  horizontal  cross  section  of  subgrid  TKE 
(Fig.  12)  shows  a  contiguous  region  of  enhanced  TKE 
with  a  relatively  large  area.  The  TKE  held  occurs  in 
coherent  bands  that  are  oriented  along  the  wave  phase 
lines  seen  in  Fig.  9.  Therefore,  the  geometry  of  the 
volumes  of  enhanced  TKE  are  complicated  in  shape, 
with  larger  horizontal  extent  than  vertical  and  an  aspect 
ratio  that  is  qualitatively  similar  to  that  of  the  inertia- 
gravity  waves. 

d.  Sensitivity  of  the  modeled  gravity  waves  to 

resolution 

To  determine  the  robustness  of  the  simulated  gravity 
waves,  numerous  resolution  sensitivity  studies  were 
completed  by  using  only  some  of  the  coarser  domains. 
For  example,  one  simulation  used  domain  1  with  no 
further  nesting  (experiment  1)  and  was  integrated  until 
0600  UTC  18  February  2001.  To  test  the  effect  of  dou- 
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(a)  Potential  Temperoture  (K)  0600  UTC  18  Feb  2001  EX-1 


0  200  400  600 

Distance  (km) 

Fig.  13.  Cross  sections  of  (a)  potential  temperature  and  (b)  wind 
speed  from  CHM  experiment  1  (domain  1)  at  0600  UTC  18  Feb  2001. 
The  cross  section  is  between  the  southwest  and  northeast  corners  of 
the  model  domain,  only  encompassing  the  horizontal  extent  of  domain 
3  (see  Fig.  6).  Also  shown  in  (a)  is  the  2-PVU  contour  of  potential 
vorticity  (thick  line). 


bling  both  the  vertical  and  horizontal  resolution,  ex¬ 
periment  2  incorporated  domains  1  and  2  until  0600 
UTC.  The  effect  of  doubling  the  vertical  resolution 
again  was  further  investigated  in  experiment  3,  which 
used  domains  1-3.  All  of  these  experiments  produced 
qualitatively  similar  results,  with  no  major  differences 
in  the  properties  of  the  jet,  upper  front,  or  gravity  waves. 

For  example,  consider  the  results  from  experiment  1 
at  0600  UTC  (Fig.  13),  through  the  same  cross  section 
as  Fig.  7  (experiment  3).  This  comparison  illustrates  the 
effect  of  doubling  the  horizontal  resolution  and  quad- 
rupoling  the  vertical  resolution.  As  expected,  experi¬ 
ment  3  shows  stronger  gradients  in  potential  tempera¬ 
ture,  stronger  (horizontal  and  vertical)  wind  shears,  and 
a  well-dehned  secondary  tropopause  fold.  Also,  the 
gravity  waves  seem  better  resolved  in  experiment  3, 
with  the  perturbations  in  potential  temperature  being 
smoother.  There  are  minor  differences  in  all  helds;  how¬ 


ever,  the  overall  dynamics  are  the  same.  Therefore,  it 
seems  apparent  that  the  results  presented  earlier  from 
experiment  3  are  relatively  insensitive  to  resolution. 

Previous  studies  (e.g.,  Lindzen  and  Fox-Rabinovitz 
1989;  Pecnick  and  Keyser  1989;  Persson  and  Warner 
1991)  have  discussed  the  importance  of  consistent  ver¬ 
tical  and  horizontal  resolution  in  the  modeling  of  at¬ 
mospheric  flows,  especially  frontogenesis.  Following 
Persson  and  Warner,  we  define  a  ratio 


AZ/A// 


where  S  is  the  slope  of  the  front  (vertical/horizontal), 
AZ  is  the  vertical  grid  spacing,  and  A//  is  the  effective 
horizontal  grid  spacing  in  the  cross-frontal  direction.  As 
shown  by  the  previous  studies,  it  is  desirable  that  the 
front  is  resolved  equally  well  in  the  horizontal  and  ver¬ 
tical,  that  is,  A  ~  1.  Otherwise,  spurious  waves  may 
result.  Persson  and  Warner  showed  that  these  spurious 
waves  form  on  scales  close  to  the  grid  scale. 

Figure  8  shows  that  in  this  case  the  upper-level  front 
is  aligned  northwest  to  southeast,  and  therefore  the  ef¬ 
fective  cross-frontal  grid  spacing  becomes  V2  X  AT. 
As  discussed  earlier  the  modeled  frontal  zone  has  a  slope 
S  ~  1/60  (vertical/horizontal),  in  the  cross-frontal  di¬ 
rection.  The  analysis  has  focused  on  domain  3,  within 
which  most  of  the  upper-level  frontogenesis  occurs.  Do¬ 
main  3  has  \H  =  \/l  X  3  km,  AZ  =  100  m,  and, 
therefore,  A  =  1.4,  which  is  close  to  its  most  desirable 
value.  Also,  given  that  the  results  from  the  CHM  are 
also  insensitive  to  resolution,  it  seems  reasonable  to 
assume  that  the  modeled  inertia-gravity  waves  pre¬ 
sented  here  are  not  spurious. 

Although  the  intensification  of  the  inertia-gravity 
waves  coincides  with  upper-level  frontogenesis  (in  do¬ 
main  3  of  the  CHM),  the  waves  may  actually  originate 
earlier  from  COAMPS.  If  that  is  the  case,  would  the 
waves  be  spurious?  CO  AMPS  has  A//  =  V2  X  1 8  km, 
AZ  =  400  m,  and,  therefore,  A  ~  0.94,  which  is  even 
closer  to  its  most  desirable  value.  Thus,  if  the  gravity 
waves  did  originate  from  CO  AMPS,  they  are  not  likely 
to  be  spurious  either. 

5.  Analysis  of  the  dropsondes  for  gravity  waves 

In  this  section,  the  NOAA  G-IV  dropsondes  are  an¬ 
alyzed  for  inertia-gravity  waves  to  determine  whether 
the  waves  in  the  simulations  are  consistent  with  the 
observations. 

The  close  spacing  of  the  dropsondes  shown  in  Fig.  2 
allows  the  detailed  mesoscale  structure  of  the  jet/front 
to  be  elucidated,  including  intense  horizontal  shear 
zones,  and  possible  coherent  phase  structures  of  gravity 
waves.  However,  the  spacing  of  the  dropsondes  is  too 
large  to  properly  sample  waves  with  horizontal  wave¬ 
lengths  of  approximately  100  km,  and  the  release  height 
of  the  sondes  is  not  high  enough  to  identify  coherent 
phase  lines  in  the  lower  stratosphere  such  as  those  in 
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Fig.  14.  Perturbation  hodographs  of  the  first  five  dropsondes  (with  good  wind  data)  released 
from  the  G-IV.  (Perturbation  velocity  is  in  m  and  the  vertical  locations  of  the  points  on  the 
hodograph  are  marked  in  km.)  Also,  the  relative  locations  of  these  dropsondes  along  the  G-IV 
flight  leg  are  described  in  Table  2. 


Fig.  7a.  Nonetheless,  the  dropsondes  do  have  high  ver¬ 
tical  resolution,  and  hodograph  techniques  (e.g..  Guest 
et  al.  2000  and  references  therein)  can  be  utilized  to 
identify  the  presence  of  an  inertia-gravity  wave  and  its 
characteristics.  However,  given  the  limited  stratospheric 
data  from  the  dropsondes,  and  the  inaccuracies  of  the 
modeled  flow  mentioned  in  the  previous  sections,  it  is 
overly  ambitious  to  compare  quantitative  details  of  the 
observed  and  modeled  inertia-gravity  waves.  Never¬ 
theless,  a  qualitative  comparison  between  the  model  and 
the  derived  wave  parameters  is  useful  for  assessing  the 
model’s  performance  in  producing  realistic  results. 

The  13  dropsondes  with  complete  wind  measure¬ 
ments  in  the  lower  stratosphere  are  analyzed  for  gravity 
waves  using  the  technique  described  by  Plougonven  et 
al.  (2003).  The  vertical  profiles  of  horizontal  velocity 
from  the  dropsondes  are  separated  into  slowly  varying 
background  profiles  and  perturbation  profiles  (that  con- 

Table  2.  Summary  of  wave  characteristics  derived  from  the  ho¬ 
dograph  method:  sonde  number,  release  time  (UTC  17  Feb  2001), 
horizontal  location  of  profile  in  Fig.  2b  (km),  vertical  wavelength 
(V),  and  horizontal  wavelength  (X^). 


Sonde 

Time 

Location 

X,  (km) 

(km) 

1 

2326 

0 

1.4-1. 8 

70-110 

2 

2331 

55 

1. 3-1.5 

105-120 

3 

2335 

110 

1.0-1. 4 

110-235 

4 

2345 

220 

1. 2-2.1 

175-310 

5 

2348 

260 

2.2-23 

105-110 

tain  vertical  scales  between  500  m  and  3.5  km)  using 
the  Alter  described  by  Scavuzzo  et  al.  (1998).  These 
perturbation  velocity  profiles  are  used  to  construct  per¬ 
turbation  hodographs.  As  described  in  previous  studies 
(e.g..  Guest  et  al.  2000;  Plougonven  et  al.  2003),  linear 
theory  of  gravity  waves  in  a  uniform  background  flow 
predicts  that  one  vertical  wavelength  of  an  inertia-grav¬ 
ity  wave  traces  an  elliptical  (perturbation)  hodograph, 
which  rotates  anticyclonically  for  an  upward-propagat¬ 
ing  wave  and  cyclonically  for  a  downward-propagating 
wave.  The  orientation  of  the  major  axis  of  the  ellipse 
is  aligned  with  the  horizontal  wave  vector  of  the  inertia- 
gravity  wave,  and  the  ratio  of  the  major  and  minor  axes 
of  the  ellipse  is  proportional  to  the  wave  frequency. 
Thus,  as  the  vertical  wavelength  can  be  readily  obtained 
from  the  sounding,  the  dispersion  relation  can  be  used 
to  determine  the  remaining  wave  characteristics.  This 
technique  has  limitations  in  conditions  with  strong  back¬ 
ground  wind  shear  and  is  only  appropriate  for  a  part  of 
the  wave  spectrum  (the  low-frequency,  long-wave 
limit).  See  Plougonven  et  al.  (2003)  for  the  specifics  of 
the  technique  and  details  of  the  limitations. 

The  hodograph  technique  is  applied  above  8  km  to 
those  dropsondes  with  good  wind  data  (above  8  km) 
shown  in  Fig.  2b.  Of  these  13  sondes,  only  the  first  5 
show  clear  evidence  of  anticyclonically  rotating  hodo¬ 
graphs.  These  hodographs  are  shown  in  Fig.  14,  and  the 
characteristics  inferred  from  these  hodographs  are  sum¬ 
marized  in  Table  2.  As  can  be  seen  from  Table  2,  these 
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five  sondes  provide  estimates  of  horizontal  and  vertical 
wavelengths  that  are  similar  to  those  determined  from 
the  numerical  simulations.  However,  the  estimates  of 
vertical  wavelengths  from  the  sondes  are  generally  be¬ 
tween  0  and  1.5  km  smaller  than  those  in  the  model. 
The  range  of  horizontal  wavelengths  surround  the  mod¬ 
eled  wavelength  of  180  km.  However,  the  limited  depth 
of  the  stratosphere  observed  by  the  dropsondes,  means 
that  only  one  or  two  vertical  wavelengths  are  sampled, 
which  adds  some  uncertainty  to  the  analysis. 

Figure  14  shows  that  between  sondes  1  and  5,  the 
orientation  of  the  wave  vector  undergoes  a  transition 
from  pointing  to  the  southwest  (or  northeast)  (sonde  1) 
to  pointing  to  the  northwest  (or  southeast)  (sonde  5).  In 
the  model,  the  wave  vector  points  strictly  to  the  south¬ 
west.  Sondes  1  and  2  show  a  wave  vector  orientation 
that  is  consistent  with  the  model,  whereas  sondes  3-5 
show  an  inconsistent  orientation.  As  seen  in  Fig.  2b, 
this  change  in  orientation  in  the  wave  vector  also  co¬ 
incides  with  the  sondes  that  passed  through  more  intense 
wind  shear.  There  are  at  least  two  important  difficulties 
in  applying  the  hodograph  method  in  regions  of  strong 
shear.  First,  Hines  (1989)  argues  that  the  orientation  of 
the  hodograph-derived  wave  vector  can  be  skewed  to¬ 
ward  the  shear  vector  by  vertical  advection  of  the  back¬ 
ground  wind  shear.  Indeed,  the  wave  vector  orientation 
of  sonde  5  is  aligned  with  the  background  shear.  Second, 
in  regions  of  stronger  shear,  the  filtering  applied  to  the 
wind  profiles  is  less  effective.  Subsequently,  signals  of 
the  strong  shear  can  be  unintentionally  incorporated  into 
the  perturbation  hodograph,  aligning  the  major  axis  of 
the  hodograph  with  the  shear  vector.  Therefore,  because 
of  these  limitations,  the  apparent  inconsistencies  in  the 
wave  vector  orientations  may  be  due  solely  to  the  lim¬ 
itations  of  the  technique  in  regions  of  strong  shear. 

The  horizontal  wavelengths  determined  using  the  ho¬ 
dograph  method  and  those  inferred  from  the  CHM  sim¬ 
ulation  show  good  agreement.  Also,  those  sondes  with 
strong  wave  signatures  are  in  the  same  region  relative 
to  the  jet  as  the  waves  in  the  model.  However,  given 
the  limitations  of  the  hodograph  technique  in  its  appli¬ 
cation  to  this  case,  it  is  difficult  to  make  definitive  con¬ 
clusions  about  the  consistency  between  the  observed 
and  modeled  waves.  Certainly,  there  appear  to  be  no 
significant  inconsistencies. 

6.  Summary  and  conclusions 

An  intense  jet  stream/upper-level  front  was  observed 
on  18  February  2001  during  SCATCAT.  These  obser¬ 
vations  were  17  closely  spaced  dropsondes  (~40  km 
separation)  released  from  NOAA’s  G-IV  weather  recon¬ 
naissance  aircraft  above  the  jet/front  system.  Following 
the  dropsonde  deployment,  the  G-IV  completed  stacked 
flight  legs  in  the  strong  vertical  shear  region  above  the 
jet  core.  While  completing  these  legs,  the  G-IV  en¬ 
countered  clear-air  turbulence  in  regions  that  were  iden¬ 
tified  as  possible  regions  of  gravity  wave  activity. 


Analysis  of  the  dropsonde  observations  showed  that 
the  maximum  wind  speed  of  the  jet  exceeded  100  m  s  ‘, 
with  very  strong  wind  shear  above  and  below  the  jet 
core.  Associated  with  this  jet  was  a  strong  upper-level 
front,  and  deep  tropopause  fold,  elucidated  by  a 
COAMPS  forecast.  The  COAMPS  forecast  compared 
qualitatively  well  with  the  observations  and  was  used 
as  the  initial  conditions  and  boundary  conditions  for  a 
finer-scale  multiply  nested  model  (CHM). 

Using  nesting,  the  CHM  simulation  resolved  the  jet/ 
front  system  on  a  grid  that  had  horizontal  grid  spacing 
of  3  km  and  vertical  grid  spacing  of  100  m.  As  the 
upper  front  intensified,  inertia-gravity  waves  intensified 
also  and  propagated  vertically  through  the  strong  neg¬ 
ative  (speed)  shear  region  above  and  to  the  anticyclonic 
side  of  the  jet  core.  These  waves  had  a  horizontal  wave¬ 
length  of  about  180  km  and  a  vertical  wavelength  of 
about  2.5  km,  and  were  shown  to  be  insensitive  to  model 
resolution. 

In  the  region  of  the  jet,  the  Richardson  number  was 
generally  low  (Ri  ~  2),  and  as  the  inertia-gravity  waves 
propagated  through  the  strong  shear  zone,  the  wave  per¬ 
turbations  significantly  modified  the  local  Richardson 
number.  Subsequently,  along  the  phase  lines  of  the  grav¬ 
ity  waves  there  appeared  bands  of  increased  and  reduced 
(Ri  <  1)  Richardson  numbers.  These  bands  of  Ri  <  1 
were  regions  of  possible  Kelvin-Helmholtz  instability. 
Further  nesting  (with  horizontal  and  vertical  grid  spac- 
ings  of  1  km  and  50  m,  respectively)  showed  that  lo¬ 
calized  volumes  of  enhanced  resolved  and  subgrid  tur¬ 
bulent  kinetic  energy  formed  in  the  bands  of  reduced 
Richardson  numbers.  This  enhanced  turbulence  was  due 
to  the  inertia-gravity  waves  and  it  is  proposed  that  such 
a  mechanism  contributed  to  the  turbulence  encountered 
by  the  G-IV. 

Even  in  the  highest-resolution  domain,  there  was  no 
convective  overturning  on  the  scale  of  the  inertia-grav¬ 
ity  wave,  and  therefore  the  waves  were  able  to  propagate 
through  the  region  of  Kelvin-Helmholtz  instability 
without  breaking  per  se.  However,  from  this  study  it  is 
difficult  to  determine  what  effect  the  instability  has  on 
the  properties  of  the  propagating  gravity  wave,  once  it 
has  propagated  through  the  shear  layer.  This  will  be 
considered  in  future  research. 

The  hodograph  technique  was  used  to  analyze  the 
dropsondes  for  inertia-gravity  waves.  This  analysis  did 
show  the  presence  of  an  inertia-gravity  wave  above  the 
jet,  whose  characteristics  displayed  some  consistency 
with  those  of  the  simulated  waves.  However,  because 
of  the  limited  stratospheric  data,  and  the  imperfections 
of  the  hodograph  method  in  regions  of  strong  wind 
shear,  the  comparisons  with  the  simulation  were  limited. 

In  conclusion,  this  study  has  examined  inertia-gravity 
waves  generated  by  a  jet  stream-upper- level  frontal  sys¬ 
tem  that  propagated  through  the  strongly  sheared  flow 
above  the  jet,  and  subsequently  generated  CAT.  Nev¬ 
ertheless,  this  study  has  investigated  inertia-gravity 
wave  and  turbulence  generation  in  only  one  part  of  a 
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baroclinic  wave:  where  inertia-gravity  waves  exist  in 
the  northwesterly  jet  stream  flow  toward  a  trough,  with 
phase  lines  that  are  closely  aligned  with  the  flow.  This 
flow  regime  was  also  studied  by  O’Sullivan  and  Dunk- 
erton  (1995).  Another  regime  was  examined  by  Zhang 
(2004):  the  southwesterly  jet  stream  flow  toward  a  ridge, 
which  generated  gravity  waves  with  phase  lines  per¬ 
pendicular  to  the  mean  flow.  Idealized  studies  (e.g., 
Zhang  2004)  can  be  used  to  examine  the  difference 
between  these  flow  regimes,  the  mechanisms  controlling 
the  wave  source,  and  the  details  of  the  wave  propaga¬ 
tion;  such  idealized  studies  are  clearer  and  less  com¬ 
plicated  than  real  cases,  such  as  the  one  considered  here. 
Idealized  studies  of  these  different  flow  regimes  and 
their  relationship  to  clear-air  turbulence  is  a  topic  of 
continuing  research. 
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